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LONG-TERM  GOALS 

The  goal  of  this  research  is  to  understand  the  effects  of  range-dependent  sediment  properties  on  the 
acoustic  field  in  2-D  shallow  water  environments.  This  information,  in  part,  is  required  to  solve  the 
statistical  inference  problem  in  inhomogeneous  shallow  water  environments. 

OBJECTIVES 

The  spatio-temporal  variability  of  the  water  column  needs  to  be  properly  accounted  for  while  looking 
for  the  effect  of  range-dependent  sediment  properties  on  the  acoustic  field.  The  objective  of  the  current 
work  is  to  create  a  high-fidelity  model  of  the  water  column  along  a  single  propagation  track  for  the 
Shallow  Water  2006  (SW06)  experiment.  The  water  column  model  is  constructed  from  oceanographic 
data  measured  at  moorings  located  along  the  acoustic  propagation  track.  The  water  column  model 
describes  oceanographic  features  on  several  space-time  scales,  including  those  ranging  from  frontal 
boundaries  to  nonlinear  internal  waves.  The  purpose  of  the  water  column  model  is  to  reduce  the 
mismatch  between  measured  and  modeled  acoustic  data  so  that  other  propagation  effects,  such  as 
range-dependent  sediment  properties  can  be  investigated  in  the  future. 

APPROACH 

The  main  goal  of  the  work  is  to  investigate  the  impact  of  range-dependent  sediment  properties  on  the 
acoustic  field  in  the  SW06  environment.  However,  during  an  analysis  of  acoustic  data  in  the  region  of 
interest,  it  was  observed  that  the  passage  of  internal  waves  caused  significant  temporal  fluctuations  in 
the  acoustic  data.  The  temporal  fluctuations  on  a  single  hydrophone  ranged  anywhere  from  a  few  dB  to 
more  than  25  dB.  This  observation  is  common  to  many  shallow  water  regions  of  the  ocean. 1’2,3,4’5 
When  trying  to  extract  information  about  range-dependent  sediment  properties  from  acoustic  data  in 
shallow  water  waveguides,  it  is  important  to  separate,  insofar  as  possible,  the  part  of  the  acoustic  field 
that  depends  on  spatio-temporal  fluctuations  of  the  water  column  from  the  part  of  the  acoustic  field  that 
depends  on  range-dependent  sediment  properties. 
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The  approach  taken  here  was  to  develop  a  data-driven  water  column  model  of  the  SW06  region  and  link 
it  to  an  acoustic  propagation  model.  The  acoustic  predictions  from  the  model  are  then  compared  with 
measured  acoustic  data  from  the  SW06  experiment  to  determine  the  utility  of  the  water  column  model, 
and  the  extent  to  which  it  can  be  used  to  predict  fluctuations  in  received  acoustic  intensity.  The  water 
column  model  encompasses  a  30  km  track  between  SW47  and  SW54  as  shown  in  Fig.  1. 


05)  •  Thermistor  a  Hydrophone  ★  Acoustic  source 


Figure  1:  (a)  Map  of  the  SW06  experimental  region  with  overlaid  bathymetric  contours  and  mooring 
locations  shown  by  dots.  The  focus  of  this  work  is  a  region  defined  by  the  line  connecting  SW47  and 
SW54.  (b)  Location  of  thermistors,  hydrophones,  and  the  acoustic  source  in  the  range-depth  plane 

between  SW47  and  SW54. 

A  224  Hz  acoustic  source  was  moored  at  SW47  and  acoustic  intensity  was  received  on  a  vertical  line 
array  (VLA)  of  hydrophones  at  SW54.  The  received  acoustic  intensity  as  a  function  of  time  was 
integrated  over  the  vertical  aperture  to  create  depth-integrated  intensity  Iz(t).  Because  Iz(t)  captures  the 
vertical  distribution  of  energy  in  the  water  column,  fluctuations  in  this  quantity  reflect  changes  in  the  net 
attenuation  experienced  by  the  acoustic  wave  as  it  undergoes  modal  scattering  by  the  internal  wave 
field.  Depth-integrated  intensity  was  chosen  as  the  acoustic  metric  for  the  comparison  between 
measured  and  modeled  data.  The  usefulness  of  the  water  column  model  will  be  assessed  by  its  ability  to 
predict  the  acoustic  depth-integrated  intensity  received  on  SW54. 

The  key  to  creating  a  useful  water  column  sound  speed  profile  from  thermistor  data  measured  at 
discrete  spatial  locations  is  the  method  by  which  the  data  are  processed  and  interpolated  in  both  space 
and  time.  Figure  2  shows  the  processing  steps  employed  to  create  a  spatio-temporal  water  column 
sound  speed  profile  from  salinity  and  temperature  recordings  at  multiple  discrete  spatial  positions.  This 
process  is  discussed  in  detail  in  Chs.  3  and  4  of  Ref.  6.  One  important  feature  of  the  process  shown  in 
Fig.  2  and  discussed  in  Ref.  6  is  that  any  number  of  oceanographic  moorings  can  be  incorporated  into 
the  oceanographic  model.  If  more  than  one  oceanographic  mooring  is  included  in  the  processing,  then 
the  internal  wave  packet  will  evolve  as  a  function  of  space  and  time.  For  this  reason,  this  method  was 
called  the  evolutionary  propagated  thermistor  string  (EPTS)  method.  If  only  a  single  oceanographic 
mooring  is  included,  then  the  oceanographic  model  resembles  the  propagated  thermistor  string  (PTS) 
method  used  in  prior  work.7,8 

The  modeled  water  column  sound  speed  profile  c(r,z,t)  was  linked  to  an  acoustic  propagation  model  to 
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Sym.  Quantity  Units 

S  salinity  [ppt] 

T  temperature  [°C] 
c  sound  speed  [m/s] 
A]  I W mode  1  amp.  [m] 
<D,  IW  mode  1  shape 
>1  IW  displacement  [m] 
i  SW47,  SW29, 

SW03,  SW54 

—  the  overbar  means  that 
a  24-hr  moving-average 
filter  has  been  applied 


Figure  2:  Data  flow  diagram  showing  how  temperature  and  salinity  data  from  discrete  spatial 
locations  are  processed  to  produce  a  water  column  sound  speed  profile. 

predict  the  acoustic  field  at  the  SW54  VLA.  The  range-dependent  acoustic  model  (RAM)  was  selected 
for  this  task.  Benchmarking  against  a  more  accurate  two-way  coupled  mode  algorithm  showed  that  the 
RAM  model  possessed  sufficient  accuracy  to  compute  the  acoustic  field  in  the  presence  of  nonlinear 
internal  waves.  This  is  important  because  the  computational  efficiency  of  RAM  made  it  possible  to 
predict  the  acoustic  field  every  15  seconds  for  more  than  3  weeks  of  the  SW06  experiment. 

Finally,  measured  and  modeled  Iz{t)  were  compared  at  the  SW54  VLA  to  assess  the  usefulness  of  the 
oceanographic  model  for  predicting  acoustic  intensity  fluctuations. 

WORK  COMPLETED 

The  work  completed  includes  the  following  items: 

1.  Creating  the  EPTS  method  for  computing  c(r,z,t )  from  measured  oceanographic  data  at  discrete 
mooring  positions. 

2.  Linking  the  water  column  model  to  RAM,  and  computing  depth-integrated  acoustic  intensity  at 
the  SW54  VLA. 

3.  Processing  acoustic  data  at  the  SW54  VLA  over  21  days  of  the  experiment. 

4.  Comparing  the  modeled  acoustic  data  to  the  measured  acoustic  data  to  validate  the  water  column 
model. 

5.  Investigating  correlations  between  the  amplitude  of  the  internal  wave  field  and  the  statistics  of 
acoustic  intensity  in  both  measured  and  modeled  data. 
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(a)  EPTS  model  using  the  SW47,  SW29,  SW03,  and  SW54  moorings 


(b)  EPTS  model  using  only  the  SW54  mooring 
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Figure  3:  Snapshots  shown  every  2  hrs  ofc(r,z,t)  generated  by  the  EPTS  model  for  internal  wave 
event  6  using  (a)  four  and  (b)  one  mooring  in  the  EPTS  model. 


RESULTS 

Snapshots  of  c(r,z,t ),  shown  every  2  hrs  during  internal  wave  event  6  (17  August  22:00  through 
18  August  10:00),  are  shown  in  Fig.  3.  Data  from  four  oceanographic  moorings  (SW47,  SW29,  SW03, 
and  SW54)  were  used  to  create  the  profiles  shown  Fig.  3(a),  while  data  from  a  single  mooring  (SW54) 
were  used  to  create  the  profiles  shown  in  Fig.  3(b).  Figure  3(a)  contains  a  range-dependent  background 
sound  speed  profile  c(r,z,t )  and  shows  the  evolution  of  the  internal  wave  packet  as  it  propagates  up  the 
continental  shelf.  In  contrast,  Fig.  3(b)  possesses  a  range-independent  c(r,z,t )  and  the  internal  wave 
packet  translates  but  does  not  evolve  as  it  propagates. 

While  Figs.  3(a)  and  (b)  appear  to  have  much  in  common,  the  differences  in  c(r,z,t)  can  dramatically 
affect  the  ability  to  accurately  predict  the  acoustic  field.  This  is  demonstrated  in  Fig.  4(a)  through  (c) 
which  show  Iz(t)  from  measured  data,  modeled  EPTS  data  using  four  moorings,  and  modeled  EPTS 
data  using  one  mooring,  respectively. 

In  general,  Iz{t)  from  Fig.  4(c)  is  not  as  accurate  as  f(l)  from  Fig.  4(b)  when  predicting  the  measured 
data.  The  model  errors  come  in  at  least  two  types.  The  first  type  of  error  is  visible  between  24  to  29 
August,  where  Iz{t)  in  Fig.  4(c)  is  5  to  10  dB  below  the  measured  value.  This  particular  error  results 
from  the  construction  of  a  range-independent  c{z,t)  from  only  ?sw54(z,0  and  Ssw54 {z,  t).  Between 
approximately  24  to  29  August  the  shelf-break  front  moved  inshore  of  the  SW54  mooring  bringing 
warmer,  saltier  water  near  SW54.  This  oceanographic  front  was  physically  localized  near  SW54  and 
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Figure  4:  Comparison  between  (a)  measured,  (b)  EPTS  modeled  (four  moorings ),  and  (c)  EPTS 
modeled  (one  mooring )  Iz(t )  with  internal  wave  event  numbers  shown  at  the  top. 


did  not  represent  the  background  water  column  conditions  over  the  entire  30  km  track.  However,  the 
sound  speed  reconstruction  used  to  generate  the  acoustic  data  shown  in  Fig.  4(c)  has  no  information 
input  about  this  range -dependent  feature  and  reconstructs  the  entire  water  column  sound  speed  profile 
from  local  conditions.  The  net  result  is  acoustic  predictions  with  large  errors  in  the  time-averaged  value 

of  hit). 

The  second  type  of  error  in  I  ft)  for  Fig.  4(c)  occurs  in  predictions  for  shorter  time  scales.  The  details  of 
the  short-time  fluctuations  in  Fig.  4(c)  do  not  reflect  those  appearing  in  Fig.  4(a).  This  type  of  error 
results  from  the  construction  of  the  internal  wave  field  from  a  single  thermistor  string.  Thus,  the  details 
of  the  internal  wave  field  are  not  as  accurate  as  they  are  for  the  EPTS  model  when  using  multiple 
moorings.  For  example,  an  inspection  of  events  12  to  13  and  18  to  19  in  Fig.  4  shows  that  the 
single-mooring  model  does  not  capture  the  approximately  10  dB  fluctuations  present  in  both  Fig.  4(a) 
and  (b).  Another  example  of  this  type  of  error  is  visible  for  event  16  for  which  the  single-mooring 
model  produces  fluctuations  approximately  10  dB  too  large.  The  change  in  the  structure  of  the  internal 
wave  field  as  it  propagates  shoreward  appears  to  significantly  affect  the  accuracy  of  the  acoustic 
prediction  on  short  time  scales. 

It  was  shown  that  large  fluctuations  in  acoustic  intensity  can  be  accurately  predicted  by  the  EPTS  model 
when  data  from  multiple  moorings  are  included.  The  ability  to  improve  the  comparison  between 
measured  and  modeled  acoustic  data  is  important  for  solving  the  range-dependent  geoacoustic  inverse 
problem.  If  the  water  column  is  not  properly  accounted  for,  then  errors  in  the  acoustic  prediction  arising 
from  errors  in  the  modeled  c(r,z,t)  will  be  translated  to  the  geoacoustic  inverse  or  inference  problem  in 
the  form  of  model  error. 
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IMPACT/APPLICATIONS 


The  impact  of  this  research  is  to  demonstrate  how  accurately  the  acoustic  field  can  be  predicted  in  the 
presence  of  internal  waves  for  long  propagation  distances  in  shallow  water.  Applications  of  this 
research  range  from  sonar  prediction  problems  to  accounting  for  a  spatio-temporally  varying  water 
column  in  inversion  and  statistical  inference  problems. 

TRANSITIONS 

The  primary  transition  for  this  project  is  a  methodology  for  reconstructing  a  modeled  sound  speed 
profile  c(r,z,t)  from  multiple  oceanographic  moorings  in  a  continental  shelf  environment. 

RELATED  PROJECTS 

None. 
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